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Abstract: We explore the sensitivity of a neutrino detector employing strongly directional high gain 
radio antennae to detect the conversion of neutrinos above 10 eV in a mountain or the earth crust. The 
directionality of the antennae will allow both, the low threshold and the suppression of background. This 
technology would have the advantage that it does not require a suitable atmosphere as optical detectors 
do and could therefore be deployed at any promising place on the planet. In particular one could choose 
suitable topographies at latitudes that are matched to promising source candidates. 
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Introduction 

Identifying the sources of the highest energy Cos- 
mic Rays (CR) is a fundamental and unsolved 
problem in astroparticle physics. If they are accel- 
erated in astrophysical objects, high energy neu- 
trino emission should be associated with this ac- 
celeration [1]. Since these neutrinos would neither 
interact with intergalactic or interstellar media nor 
be deflected by magnetic fields, neutrino astron- 
omy carries the hope of identifying these elusive 
sources of highest energy CRs. 

Various detectors have been proposed for neutrino 
energies above 10 16 eV, most of them relying on 
mountains to provide suitable target mass for v T 
with the ensuing r lepton decay providing a de- 
tectable air shower (AS) after the r leaves the tar- 
get [2]. As high energy neutrinos get absorbed in 
the earth, geometries for such events are restricted 
to within a few degrees of horizontal [3]. 

Optical detectors exploiting the fluorescence emis- 
sion of horizontal showers as well scintillator tele- 
scopes have been proposed to measure the promis- 
ing v T signature [4] . Here we explore what could 
be a simpler and cheaper technology in the field: 
strongly directional (high gain, low threshold) ra- 
dio antennas. All they require is a rigid mechani- 
cal support structure and like all the other detectors 



some electronics. But no optical components have 
to be calibrated and maintained, and the detector 
could be made largely insensitive to the local cli- 
mate at a promising detector site. The challenge 
for our detector will lie in isolating the signal. 

Geosynchrotron radio signals 

AS produce geosynchrotron radio emission due to 
the deflection of secondary shower electrons and 
positrons in the earth's magnetic field [5]. In this 
study for the first time such radio signals are stud- 
ied from horizontal showers. Using the REAS2 
Monte Carlo code [6] we calculate the radio emis- 
sion from 10 15 eV and 10 16 eV horizontal, 7r + 
induced AS propagating through a constant den- 
sity atmosphere of 1.058 x 1CP 3 g cm~ 3 , cor- 
responding to approximately 1500 m above sea 
level. The magnetic field was set to a strength 
of 0.4 Gauss with an inclination of 60°, a conser- 
vative estimate for Central Europe and Northern 
America. The AS were simulated as propagating 
from north to south; the corresponding geomag- 
netic angle is thus 60°. This is the most conser- 
vative choice, as it effectively minimizes the field 
component perpendicular to the propagation direc- 
tion of the AS. Atmospheric curvature can be ne- 
glected here, as the AS need only a few kilometers 
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to develop in the lower atmosphere. The AS sim- 
ulations were carried out using CORSIKA 6.502 
[7], with shower-to-shower fluctuations taken out 
by simulating 25 air showers per parameter set and 
then selecting a shower with a typical longitudinal 
evolution profile for the radio simulations. 
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Figure 1: Lateral dependence of geosynchrotron 
radio emission field strenghts of 10 16 eV hori- 
zontal air showers for different observer distances 
from the r decay point. 10 15 eV values are lower 
by approximately a factor of 10. 

One important question is which observer distance 
is suited best for radio observations of these hori- 
zontal AS. Figure 1 illustrates how the lateral pro- 
file of the radio signal spreads out as the distance 
from the starting point of the AS increases. Despite 
strong relativistic beaming of geosynchrotron radi- 
ation, the lateral slope changes significantly with 
the distance of the observer from the r decay point. 
An observing distance of 20 km seems to consti- 
tute a good compromise: the signal in the center 
region is still strong while its lateral spread does 
not require too dense an array of antennae. 

Also important is the selection of a suitable obser- 
vation frequency band. This clearly depends on 
the actual noise situation at the selected observ- 
ing site. In the absence of man-made radio in- 
terference, atmospheric and galactic noise set the 
limits. In Fig. 2 we show what signal-to-noise ra- 
tios (SNR, defined as peak power of the signal di- 
vided by power of the noise in the band of inter- 
est) can be expected for 10 16 eV AS at 20 km ob- 
serving distance and 225 m vertical observer off- 
set from the shower axis. The values are calcu- 
lated for an isotropic radiator and a combination 
of day-time atmospheric noise and galactic noise 
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Figure 2: Expected SNR as a function of observ- 
ing frequency band for 10 16 eV showers initiated at 
20 km distance and 225 m vertical observer offset 
from the shower axis, calculated for an isotropic ra- 
diator. 10 15 eV values are lower by approximately 
a factor of 100. 

based on measurements by the CCIR/ITU-R (Inter- 
national Telecommunication Union). As geosyn- 
chrotron emission has a steeply falling frequency 
spectrum, it is important to include low frequen- 
cies. As can be seen in Fig. 2 an observing band- 
width from 20 to 80 MHz would be desirable. Be- 
low 20 MHz, atmospheric noise gets very strong, 
and above 80 MHz FM radio transmitters could 
pose problems. 
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Figure 3: Expected SNR as a function of observ- 
ing frequency band for 10 16 eV showers at 20 km 
distance, calculated for an antenna with 16 dBi an- 
tenna gain. 10 15 eV values are lower by approxi- 
mately a factor of 100. 

SNR of order unity are too low for self-triggered 
measurements of geosynchroton radiation. Using 
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directional antennae pointing at the target mass 
will significantly improve the SNR. For the en- 
visioned broad-band measurements, logarithmic- 
periodic dipole antennae (LPDAs) can achieve an- 
tenna gains of 10 dBi. If three such antennae are 
phase coupled, effective gains of up to 16 dBi can 
be reached. This boosts the SNR into a region 
where measurements seem feasible up to axis off- 
sets of ~ 300 m, as illustrated in Fig. 3. 

For the specific relative geometry of AS and mag- 
netic field explored here the lateral distribution of 
the radio signal is very different along the hori- 
zonatal and the vertical axis. As multiple mea- 
surements along the vertical also allow to place 
additional constraints on the zenith angle of the 
observed shower and as the signal falls off more 
slowly along the vertical axis, more than one an- 
tenna should be used along the vertical direction. 

The SNR required for self-triggered operation of 
an antenna array depends on many factors such as 
the multiplicity of antennae used in coincidence 
and the required detection efficiencies. These is- 
sues would have to be addressed in a specific pro- 
posal for such a detector. 

Neutrino event rates 

A threshold close to 10 16 eV for horizontal AS 
is well matched to the effective threshold beyond 
which r leptons can escape from significant depth 
inside a rock mass: at 10 15 eV the decay length of 
a r (in vacuum) is only ~50 m. At 10 17 eV this de- 
cay length will have grown to 5 km, which allows 
for r from a reasonable amount of target mass to 
escape the rock before they decay, and on the other 
hand constrains the decay volume needed behind 
the rock to be reasonably small. 

Using a modified ANIS [8] code we examine the 
probability that a v T interacts in rock with the ensu- 
ing t decay initiating an AS outside of the rock. In 
ANIS we use the option of CTEQ5 deduced cross- 
sections. A smooth approximation for the energy 
loss of the t leptons in rock and air is included in 
the calculation. Depending on the energy of the in- 
coming v T , an equilibrium is reached at some point 
inside the rock between v T interactions producing 
new r leptons and the decay of r leptons that were 
produced in neutrino interactions further upstream. 



Our simulations show that for standard rock this 
equilibrium is reached after 1 km at 10 15 eV, 2 km 
at 10 16 eV, and 10 km at 10 17 eV. 

25 M v T were injected at each of the energies 
mentioned above. From the above simulations of 
radio signals we infer an effective threshold of 
8 x 10 15 eV for AS detection to estimate how many 
neutrinos we might detect. Fig. 4 shows the dis- 
tributions of decay vertices above that energy from 
the 10 16 eV and 10 17 eV simulations. 
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Figure 4: Distribution of decay vertices that result 
in showers with an energy above 8x 10 15 eV. The 
blue histogram is for a v T energy of 10 17 eV, the 
red one for 10 16 eV. 

Within the first 10 km after the mountain we see 
less than 200 decays for the lower v T energy and 
about 40,000 for the higher one. At higher ener- 
gies efficiency will be affected by r decaying be- 
hind the detector "volume". The detection effi- 
ciency for 10 17 eV v r is e = 0.0016. A detector el- 
ement of two stations separated vertically by 200 m 
would allow to collect data over a vertical range of 
400 m. Working with a 200m horizontal displace- 
ment between such detector elements one would 
hope to have a decent efficiency for threefold coin- 
cidences over an area of roughly 160,000 m 2 for 
two such elements, and each vertical expansion 
of the array by one additional element would add 
roughly 80,000 m 2 . The exposure of an eight an- 
tenna array accumulated over one year would be 
31,536,000 s x 320,000 m 2 w 10 17 cm 2 s for 
a pointlike source. As the LDPAs typically have 
an opening angle of 65°, detectors situated 30 km 
from a 1 .5 km high mountain range would cover a 
solid angle of 5.7 x 10 -2 sr. 



Radio detection of neutrinos 



A source would only be seen if it is behind a tar- 
get mass at the local horizon of the detector. The 
time any given source spends near the horizon de- 
pends on its declination. Figure 5 summarizes the 
total observation times that can be expected if the 
source can be followed within ±1.5° of the hori- 
zon and gives the azimuthal coverage required to 
follow the source along that segment of its path. If 
a source were to just rise or set vertically through a 
±1.5° detector aperture, it would be visible for less 
than 1% of the total time. If the location was cho- 
sen to accommodate the requisite azimuthal cover- 
age shown in Fig. 5, then one would expect to "see" 
the source about 15% of total time for a latitude of 
45°. The proposed antennae will cover the requi- 
site 30° in azimuth. Putting all this together, if we 
can measure v T between 10 16 eV and 10 17 eV, with 
the efficiencies estimated above and linearly inter- 
polated between the two energies, it would take an 
array of a little more than 130k antennae to observe 
one event per year from Galactic sources [9]. 





Figure 5: Upper panel: Percent of total time spent 
in a ±1.5° band around the horizon if the upper 
or lower culmination point is arranged to be 1.5° 
above or below the horizon. Lower panel: total 
number of degrees that have to be covered in az- 
imuth in order to attain that maximal time. 

To give an example of the sensitivity reach to 
the diffuse neutrino flux we consider, <f> Vr ~ 
1CT 3 (E u /GcV)- 2 - 54 GcV" 1 cm" 2 s" 1 sr" 1 , 
which is expected if extragalactic cosmic rays 
(from transparent sources) begin dominating 
the observed spectrum at energies as low as 



~ 10 17 - 6 eV [10], as suggested by recent HiRes 
data [11]. For such a ^ T -flux, the expected number 
of events (with E v > 10 17 eV) per year per eight 
antennae is about 0.003. 



Conclusions 

The energy range at which the detector works 
is well matched to the problem of v T detection 
through r decay in the atmosphere. Cosmic ray de- 
tectors with upward of a thousand detector stations 
are proven, keeping open in principle the possibil- 
ity of deploying an antenna array big enough for 
the isotropic fluxes. To further pursue this route, 
emphasis will have to be put on developing inex- 
pensive and reliable detector stations for such a de- 
tector. Radio detection carries that promise. 
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